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19  ABSTRACT 

20  The  dataset  from  the  SUPERTANK  laboratory  experiment  was  analyzed  to 

21  examine  wave  runup  and  the  corresponding  upper  limit  of  beach-profile  change.  Thirty 

22  SUPERTANK  runs  were  investigated  that  included  both  erosional  and  accretionary  wave 

23  conditions  with  random  and  monochromatic  waves.  The  upper  limit  of  beach  change  Ul 

24  was  found  to  approximately  equal  the  vertical  excursion  of  total  wave  runup,  R^.  An 

25  exception  was  runs  where  beach  or  dune  scarps  were  produced,  which  substantially  limit 

26  the  uprush  of  swash  motion  to  produce  a  much  reduced  total  runup.  Based  on  the 

27  SUPERTANK  dataset,  the  vertical  extent  of  wave  runup  above  mean  water  level  on  a 

28  beach  without  scarp  fonnation  was  found  to  approximately  equal  the  significant  breaking 

29  wave  height,  //bs.  Therefore,  a  new  and  simple  relation  Rtw  =  Hhs  is  proposed.  The 

30  linear  relationship  between  total  runup  and  breaking  wave  height  is  supported  by  a 

3 1  conceptual  derivation.  In  addition,  the  relation  is  extended  to  UL  =  Rtw  =  Hhs ,  to 

32  approximate  the  upper  limit  of  beach  change.  This  formula  accurately  reproduced  the 

33  measured  upper  limit  of  beach  change  from  the  three-dimensional  experiments  in  the 

34  Corps’  Large-scale  Sediment  Transport  Facility.  For  the  studied  laboratory  cases, 

35  predictions  of  wave  runup  were  not  improved  by  including  a  slope-dependent  surf- 

36  similarity  parameter.  The  limit  of  wave  runup  was  substantially  less  for  monochromatic 

37  waves  than  for  random  waves,  attributed  to  absence  of  low-frequency  motion. 

38  Additional  Index  Words  Beach  erosion,  nearshore  sediment  transport,  wave 

39  breaking,  cross-shore  sediment  transport,  physical  modeling,  surf  zone  processes. 

40 
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41  INTRODUCTION 

42 

43  Accurate  prediction  of  the  upper  limit  of  beach  change  is  necessary  for  assessing 

44  large  morphologic  changes  induced  by  extreme  storms.  The  upper  limit  of  beach  change 

45  is  controlled  by  wave  breaking  and  the  subsequent  wave  runup.  During  storms,  wave 

46  runup  is  superimposed  on  the  elevated  water  level  due  to  storm  surge.  Wang  et  al.  (2006) 

47  found  the  highest  elevation  of  beach  erosion  induced  by  Hurricane  Ivan  in  2004  to  be 

48  considerably  greater  than  the  measured  storm-surge  level,  indicating  that  wave  runup 

49  played  a  significant  role  in  the  upper  limit  of  beach  erosion.  The  limit  of  wave  runup  is 

50  also  a  key  parameter  in  the  application  of  the  storm-impact  scale  by  Sallenger  (2000). 

5 1  The  Sallenger  scale  categorizes  four  levels  of  morphologic  impact  by  storms  through 

52  comparison  of  the  highest  elevation  reached  by  stonn  water  (combined  stonn  surge  and 

53  wave  runup)  and  a  representative  elevation  of  the  barrier  island  (e.g.,  the  top  of  the 

54  foredune  ridge).  Quantification  of  wave  runup  and  its  relationship  to  the  upper  limit  of 

55  morphologic  change  are  required  for  understanding  and  predicting  beach-profile  changes. 

56  Wave  runup  is  composed  of  wave  setup  and  swash  runup,  defined  as  a  super- 

57  elevation  of  the  mean  water  level  and  fluctuation  about  that  mean,  respectively  (Guza 

58  and  Thornton,  1981;  Holland  et  al.,  1995;  Holman  and  Sallenger,  1985;  Nielsen, 

59  1988;  Yamamoto,  Tanimoto,  and  Harshinie,  1994).  Several  formulas  have  been 

60  developed  to  predict  wave  setup  and  runup.  Based  on  laboratory  experiments,  Hunt 

61  (1959)  proposed  various  formulas  estimating  wave  uprush,  R,  on  seawalls  and 

62  breakwaters,  and  the  “Hunt  formula”  continues  in  use: 

A  =  i  Q  tan/? 

H  '  Jh7T0 


63 


(1) 
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64  where  tan  ft  =  beach  slope,  H  =  wave  height,  typically  taken  to  be  the  deep-water  wave 

65  height,  and  L0  =  deepwater  wavelength.  Hedges  and  Mase  (2004)  modified  Hunt’s 

66  original  formula  to  include  the  contribution  of  wave  setup. 

67  Bowen,  Inman,  and  Simmons  (1968)  derived  a  wave  setup  slope  on  a  unifonnly 

68  sloping  beach  for  monochromatic  waves  as: 

69  —  =  -/6— ,  K  =  (1  +  2.67y~2)-1  (2) 

ox  ox 

70 

7 1  where  h  =  still-water  depth,  I]  =  wave  setup,  x  =  cross-shore  coordinate,  and  breaker 

72  index  y  =  H/( rj  +  h).  Based  on  both  theory  and  laboratory  measurements,  Battjes 

73  (1974)  found  the  maximum  setup  under  a  monochromatic  wave  77m  to  occur  at  the  still- 

74  water  shoreline,  77  m/M>  =  0.3y,  where  /R  =  breaking  wave  height.  Taking  the  commonly 

75  used  value  of  0.78  for  y,  the  maximum  setup  yielded  from  the  Battjes  formula  is  about 

76  23%  of  the  breaking  wave  height. 

77  The  development  of  most  formulas  predicting  the  limits  for  wave  runup  has 

78  involved  comparisons  to  field  measurements.  Based  on  measurements  made  on 

79  dissipative  beaches,  Guza  and  Thornton  (1981)  suggested  that  the  setup  at  the 

80  shoreline  is  linearly  proportional  to  the  significant  deepwater  wave  height  //0: 

81  %  =  0.17H0  (3) 

82 

83  In  a  following  study,  Guza  and  Thornton  (1982)  found  the  significant  wave  runup  Rs 

84  (including  both  wave  setup  and  swash  runup)  is  also  linearly  proportional  to  the 

85  significant  deepwater  wave  height: 

Rs  =3.48  +  0.71  H0 


86 


(units  of  centimeters) 


(4) 


87 

88  Comparing  Eqs.  (3)  and  (4),  the  entire  wave  runup  is  approximately  4  times  the 
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contribution  of  wave  setup,  i.e.,  swash  runup  constitutes  a  significant  portion, 
approximately  75%,  of  the  total  elevated  water  level.  According  to  Huntley  et  al. 
(1993),  Eq.  (4)  is  the  best  choice  for  predicting  wave  runup  on  dissipative  beaches. 

Based  on  field  measurements  on  highly  dissipative  beaches,  Ruessink,  Klenhans,  and 
van  den  Beukel  (1998)  and  Ruggiero  et  al.  (2001)  also  found  linear  relationships,  but 
with  slightly  different  empirical  coefficients. 

Based  on  field  measurements,  Holman  (1986)  and  several  similar  studies 
(Holman  and  Sallenger,  1985;  Ruggiero,  Holman,  and  Beach,  2004;  Stockdon  et 
al.,  2006)  argued  that  more  accurate  predictions  for  intermediate  beaches  can  be  obtained 
by  including  the  surf  similarity  parameter,  following  the  Hunt  formula  (Eq.  1): 


(5) 


Holman  (1986)  found  a  dependence  of  the  2%  exceedence  of  runup  /??  on  the  deepwater 
significant  wave  height  and  the  (offshore)  surf  similarity  parameter: 

R2  =(0.83£  +  0.2)//0  (6) 


Stockdon  et  al.  (2006)  expanded  upon  the  Holman  (1986)  analysis  with  additional  data 
covering  a  wider  range  of  beach  slopes  and  developed  the  empirical  equation: 


1  H,L, (0.563 tan j8{2  +0.004)  2 

R2=l.l  0.35 tan /?f(//0Z0)2  +- - - - =L_ 


(7) 


108 
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109  Realizing  the  variability  of  beach  slope  in  terms  of  both  definition  and  measurement, 

110  Stockdon  et  al.  (2006)  defined  the  foreshore  beach  slope  as  the  average  slope  over  a 

111  region  of  two  times  the  standard  deviation  of  a  continuous  water  level  record. 

112  With  the  exception  of  the  original  derivation  by  Bowen,  Inman,  and  Simmons 

113  (1968),  most  predictive  fonnulas  for  wave  runup  on  a  natural  beach  have  been 

114  empirically  derived  based  on  field  measurements  over  dissipative  and  intermediate 

115  beaches.  Field  measurements  of  wave  runup  were  typically  conducted  with  video 

116  imagery  and/or  resistance  wire  generally  5  to  20  cm  above  and  parallel  to  the  beach  face. 

117  Holland  et  al.  (1995)  concluded  that  these  two  measurement  methods  are  comparable  in 

118  producing  accurate  results. 

119  Almost  all  the  aforementioned  field  studies  focused  mainly  on  the  hydrodynamics 

120  of  wave  runup,  with  little  discussion  on  the  corresponding  morphologic  response, 

121  particularly  the  upper  limit  of  beach-profile  change.  Thus,  in  contrast  to  a  considerable 

122  number  of  studies  on  wave  runup,  data  are  scarce  that  relate  the  limit  of  wave  runup  with 

123  the  resulting  beach  change.  In  other  words,  the  limit  of  beach  change  as  related  to  wave 

124  runup  has  not  been  well  documented. 

125  Data  from  the  prototype-scale  laboratory  experiments,  including  those  conducted 

126  at  SUPERTANK  (Kraus,  Smith,  and  Sollitt,  1992;  Kraus  and  Smith,  1994)  and 

127  Large-scale  Sediment  Transport  Facility  (LSTF)  (Hamilton,  et  al.,  2001;  Wang,  Smith, 

128  and  Ebersole,  2002),  are  examined  in  this  paper  to  study  the  limit  of  wave  runup  and 

129  corresponding  limit  of  beach  or  dune  erosion.  Specifically,  this  study  examines  1)  the 

130  levels  of  total  wave  runup,  including  swash  runup  and  wave  setup;  2)  time-series  of 

131  beach-profile  change  under  erosional  and  accretionary  waves;  3)  the  relationship  between 
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132  the  waves  and  beach-profile  change;  and  4)  the  accuracy  of  existing  wave  runup 

133  prediction  methods.  A  new  empirical  formula  predicting  the  limits  of  wave  runup  and 

134  that  of  beach  change  is  proposed  based  on  the  prototype-scale  laboratory  data. 

135 

136  METHODS 

137  SUPERTANK  and  LSTF  Experiments 

138  Data  from  two  movable -bed  laboratory  studies,  SUPERTANK  and  LSTF  (Fig.  1), 

139  are  examined  to  quantify  the  upper  limits  of  beach-profile  change,  wave  runup,  and  their 

140  relationship.  Both  experiments  were  designed  to  measure  sediment  transport  and 

141  morphology  change  under  varying  prototype  wave  conditions.  Dense  instrumentation  in 

142  the  laboratory  setting  allows  for  well-controlled  and  accurate  measurement  of 

143  hydrodynamic  conditions  and  morphological  change.  SUPERTANK  was  a  two- 

144  dimensional  wave  channel  with  beach  change  induced  primarily  by  cross-shore  processes, 

145  whereas  the  LSTF  was  a  three-dimensional  wave  basin  with  both  cross-shore  and 

146  longshore  sediment  transport  inducing  beach  change. 

147  SUPERTANK  was  a  multi-institutional  effort  sponsored  by  the  U.S.  Army  Corps 

148  of  Engineers  and  conducted  at  the  O.H.  Hinsdale  Wave  Research  Laboratory  at  Oregon 

149  State  University  from  July  29  to  September  20,  1991.  This  facility  is  the  largest  wave 

150  channel  in  the  United  States  that  can  contain  a  sandy  beach  through  which  experiments 

151  comparable  to  the  magnitude  of  naturally  occurring  waves  can  be  conducted  (Kraus, 

152  Smith,  and  Sollitt,  1992).  The  SUPERTANK  experiment  measured  total-channel 

153  hydrodynamics  and  sediment  transport  along  with  the  resulting  beach-profile  change. 

154  The  wave  channel  is  104  m  long,  3.7  m  wide,  and  4.6  m  deep  (the  still  water  level  was 
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155  typically  1.5  m  below  the  top  during  SUPERTANK)  with  a  constructed  sandy  beach 

■J 

156  extending  76  m  offshore  (Fig.  1  upper).  The  beach  was  composed  of  600  m  of  fine, 

157  well-sorted  quartz  sand  with  a  median  size  of  0.22  mm  and  a  fall  speed  of  3.3  cm/s.  The 

158  wave  generator  and  wave  channel  were  equipped  with  a  sensor  to  absorb  the  energy  of 

159  reflected  waves  (Kraus  and  Smith,  1994).  The  water-level  fluctuations  were  measured 

160  with  16  resistance  and  10  capacitance  gauges.  These  26  gauges,  spaced  no  more  than  3.7 

161  m  apart,  provided  high  resolution  of  wave  propagation,  especially  in  the  swash  zone. 

162  The  beach  profile  was  surveyed  following  each  20-  to  60-min  wave  run.  The 

163  initial  profile  was  constructed  based  on  the  equilibrium  beach  profile  developed  by 

1 64  Bruun  (1954)  and  Dean  ( 1 977)  as: 

165  h(x)  =  Ax213  (8) 

166  where  h  =  still-water  depth,  x  =  horizontal  distance  from  the  shoreline,  and  A  =  a  shape 

167  parameter,  which  for  SUPERTANK  corresponded  to  a  median  grain  size  of  0.30  mm. 

168  The  initial  beach  was  built  steeper  with  a  greater  4 -value  to  ensure  adequate  water  depth 

169  in  the  offshore  area  (Wang  and  Kraus,  2005).  For  efficiency,  most  SUPERTANK  cases 

170  were  initiated  with  the  final  profile  of  the  previous  run.  Approximately  350  profile 

171  surveys  were  made  by  using  an  auto-tracking,  infrared  Geodimeter  targeting  prism 

172  attached  to  a  survey  rod  mounted  on  a  carriage  pushed  by  researchers.  Three  along- 

173  channel  lines  were  surveyed.  Only  the  center  line  was  analyzed  in  this  study.  Wave- 

174  processing  procedures  are  discussed  in  Kraus  and  Smith  (1994).  To  separate  incident- 

175  band  wave  motion  from  low-frequency  motion,  a  non-recursive,  low-pass  filter  was 

176  applied.  The  period  cutoff  for  the  filter  was  set  to  twice  the  peak  period  of  the  incident 

177  waves. 
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178  The  LSTF  is  a  three-dimensional  wave  basin  located  at  the  U.S.  Army  Corps  of 

179  Engineers  Coastal  and  Hydraulic  Laboratory  in  Vicksburg,  Mississippi.  Operation 

180  procedures  are  discussed  in  Hamilton  et  al.  (2001).  The  LSTF  was  designed  to  study 

181  longshore  sediment  transport  (Wang,  et  al.,  2002;  Wang,  Smith,  and  Ebersole,  2002). 

182  The  LSTF  is  capable  of  generating  wave  conditions  comparable  to  the  naturally  occurring 

183  wave  heights  and  periods  found  along  low-energy  open  coasts  and  bays.  The  LSTF  has 

184  dimensions  of  30  m  across-shore,  50  m  longshore,  with  walls  1.4  m  high  (Fig.  1  lower). 

3 

185  The  beach  was  designed  in  a  trapezoidal  plan  shape  composed  of  approximately  150  m 

186  of  very  well  sorted  fine  quartz  sand  with  a  median  grain  size  of  0.15  mm  and  a  fall  speed 

187  of  1 .8  cm/s.  Initial  construction  of  the  beach  was  also  based  on  the  equilibrium  profile 

188  (Eq.  8).  The  beach  profile  was  surveyed  using  an  automated  bottom-tracking  profiler 

189  capable  of  resolving  bed  ripples.  The  beach  was  typically  replenished  after  3  to  9  hours 

190  of  wave  activity.  Long-crested  and  unidirectional  irregular  waves  with  a  relatively  broad 

191  spectral  shape  were  generated  at  a  10  deg  incident  angle  in  the  horizontal  section  of  the 

192  basin.  The  wave  height  and  peak  wave  period  were  measured  with  capacitance  wave 

193  gauges  sampling  at  20  Hz,  statistical  wave  properties  were  calculated  by  spectral 

194  analysis.  The  experimental  procedures  in  LSTF  are  described  in  Wang,  Smith,  and 

195  Ebersole  (2002). 

196 

197  Data  Analysis 

198  Although  the  entire  SUPERTANK  dataset  is  available,  five  cases  with  a  total  of 

199  30  wave  runs  were  selected  from  the  20  initial  cases.  The  selection  was  based  on  the 

200  particular  purpose  of  the  wave  run,  the  trend  of  net  sediment  transport,  and  the  beach 
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201  response.  Time-series  of  beach-profile  change,  cross-shore  distribution  of  wave  height, 

202  and  mean  water  level  were  analyzed.  Scarp  presence  was  also  identified.  The  upper 

203  limits  Ul  for  the  non-scarped  beach  profile  runs  were  identified  based  on  the  upper 

204  profile  convergence  point,  above  which  no  beach  change  occurred.  The  upper  limit 

205  identified  for  the  scarped  runs  was  at  the  elevation  of  the  scarp  toe.  The  other  location 

206  examined  was  the  lower  limit  of  beach  change  Tl-  The  lower  limit,  or  lower  profile 

207  convergence  point,  was  identified  at  the  depth  contour  below  which  no  change  occurred. 

208  For  the  30  SUPERTANK  wave  runs  examined,  the  water  level  and  zero-moment 

209  wave  height  were  analyzed.  From  the  cross-shore  wave  height  distribution  (or  wave- 

210  energy  decay),  the  breaker  point  H\,  was  defined  at  the  location  with  a  sharp  decrease  in 

2 1 1  wave  height  (Wang  et  al.,  2002).  The  total  wave  runup  Rtw  was  defined  by  the  location 

212  and  beach  elevation  of  the  swash  gauge  that  contained  a  value  larger  than  zero  wave 

213  height,  i.e.,  water  reached  that  particular  gauge.  The  above  procedure  did  not  involve  any 

214  statistical  analysis,  but  rather  was  determined  by  the  measurements  available  from 

215  SUPERTANK.  Hence,  there  may  be  some  differences  between  the  RXw  determined  in  this 

216  study  and  the  2%  exceedence  of  runup  (f?2%)  as  appears  in  some  predictive  equations, 

217  obtained  from  video  (e.g.,  Holland  et  al.,  1995)  and  horizontally  elevated  wires  (e.g., 

218  Guza  and  Thornton,  1982). 

219  Two  LSTF  experiments,  one  conducted  under  random  spilling  breaker  waves  and 

220  one  under  random  plunging  breaker  waves,  were  examined  in  this  study.  The  LSTF  data 

22 1  were  examined  for  the  upper  limit  of  beach  change.  The  beach-profiles  analyzed  here 

222  were  surveyed  through  the  middle  of  the  basin.  The  maximum  runup  was  not  directly 
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measured  due  to  the  lack  of  swash  gauges.  The  main  objective  of  the  LSTF  analysis  was 
to  apply  the  SUPERTANK  results  to  a  three-dimensional  beach. 


RESULTS 


Overall,  30  SUPERTANK  wave  runs  and  two  LSTF  wave  cases  were  analyzed 
(Table  1).  The  two  LSTF  cases,  under  a  spilling  and  a  plunging  breaker,  examined  the 
effect  of  the  breaker  types  on  sediment  transport  and  the  resulting  beach-profile  change. 
The  thirty  SUPERTANK  wave  runs  are  composed  of  twelve  erosional  random  wave 
runs,  three  erosional  monochromatic  wave  runs,  seven  accretionary  random  wave  runs, 
three  accretionary  monochromatic  wave  runs,  and  five  dune  erosion  random  wave  runs. 

In  Table  1,  the  first  two  numbers  in  the  wave  run  ID  “10A60ER”  indicate  the  major  data 
collection  case,  the  letter  “A”  indicates  a  particular  wave  condition,  and  the  numerals 
indicate  the  duration  of  wave  action  in  minutes.  The  notation  used  in  Table  1  and  in  all 
equations  is  listed  in  Appendix  I.  The  erosional  and  accretionary  cases  were  designed 
based  on  the  Dean  number  N, 

N  =  (9) 

rri  v  ' 

wl 

where  H\,s  =  significant  breaking  wave  height;  w  =  fall  speed  of  the  sediment,  and  T  = 
wave  period  (Kraus,  Smith,  and  Sollitt,  1992). 

Beach-Profile  Change 
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245  Analysis  of  the  time-series  of  beach-profile  change  for  the  SUPERTANK 

246  experiment  can  be  found  in  Roberts,  Wang,  and  Kraus  (2007).  The  first 

247  SUPERTANK  wave  run,  ST_10A,  was  conducted  with  a  monotonic  initial  profile  (Eq. 

248  8).  Figure  2  illustrates  four  time-series  beach-profiles  surveyed  at  initial,  60,  130,  and 

249  270  min.  The  design  wave  conditions  are  included  as  an  inset  in  the  figures,  where  Hmo 

250  is  the  zero-moment  wave  height,  Tp  is  the  peak  spectral  wave  period,  and  n  is  spectral 

25 1  peakedness.  Significant  beach-profile  change  occurred  with  substantial  shoreline 

252  recession,  along  with  the  development  of  an  offshore  bar.  Initially,  the  foreshore 

253  exhibited  a  convex  shape  while  the  end  profile  was  concave.  The  upper  limit  of  beach- 

254  profile  change  was  measured  at  0.66  m  above  mean  water  level  (MWL)  for  all  three  time 

255  segments.  An  apparent  point  of  profile  convergence  was  measured  at  the  1.35m  depth 

256  contour,  beyond  which  profile  elevation  change  cannot  be  clearly  identified. 

257  The  subsequent  wave  runs  were  conducted  over  the  final  profile  of  the  previous 

258  wave  run,  i.e.,  over  a  barred  beach.  The  beach-profile  changes  are  detectable,  but  much 

259  more  subtle  than  the  initial  run  (Fig.  2),  especially  for  the  accretionary  wave  runs  with 

260  lower  wave  heights.  Figure  3  shows  an  example  of  an  accretionary  wave  run,  ST_30A. 

261  The  upper  limit  was  detennined  to  be  at  0.3 1  m  above  MWL  (Fig.  3  lower).  One  of  the 

262  surveys  (60  min)  exhibited  some  changes  above  that  convergence  level;  however  these 

263  changes  may  be  attributable  to  survey  error.  The  offshore-profile  convergence  point  was 

264  determined  at  a  depth  contour  of  around  1  m. 

265  A  scarp  developed  in  some  of  the  erosional  wave  runs  (Fig.  4).  The  scarp  was 

266  induced  by  wave  erosion  of  the  base  of  the  dune  or  the  dry-beach,  subsequently  causing 

267  the  overlying  sediment  to  become  unstable  and  collapse.  The  resulting  beach  slope 
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268  directly  seaward  of  the  scarp  tends  to  be  steeper  than  on  a  non-scarped  beach.  The  upper 

269  limit  of  beach  change  is  apparently  at  the  top  of  the  scarp,  controlled  by  the  elevation  of 

270  the  beach  berm  or  dune,  and  does  not  necessarily  represent  the  vertical  extent  of  wave 

27 1  action.  The  upper  limit  of  beach  change  in  this  study  was  selected  at  the  base  of  the  near- 

272  vertical  scarp,  measured  at  0.31  m  above  MWL  during  ST_60A  (Fig.  4).  Therefore,  for 

273  the  scarped  case,  the  upper  limit  was  controlled  by  both  wave  action  and  gravity-driven 

274  dune  collapse.  Little  beach-profde  change  was  observed  offshore. 

275  SUPERTANK  experiments  also  included  several  runs  with  monochromatic  waves 

276  Beach-profile  change  under  monochromatic  wave  action  was  substantially  different  from 

277  those  under  the  more  realistic  random  waves  (Fig.  5).  The  monochromatic  waves  tended 

278  to  create  erratic  and  undulating  profiles.  For  STIO,  the  upper  limit  of  beach-profile 

279  change  was  estimated  at  around  0.50  m  and  varied  slightly  during  the  different  wave  runs 

280  The  erratic  profile  evolution  did  not  seem  to  approach  a  stable  equilibrium  shape,  and  it 

28 1  did  not  have  an  apparent  profile  convergence  point.  In  addition,  the  profile  shape 

282  developed  under  monochromatic  waves  does  not  represent  profiles  typically  measured  in 

283  the  field  (Wang  and  Davis,  1998).  This  implies  that  morphological  change  measured  in 

284  movable-bed  laboratory  experiments  under  monochromatic  waves  may  not  be  applicable 

285  to  a  natural  setting. 

286  Similar  analyses  were  also  conducted  for  the  data  from  the  LSTF.  The  waves 

287  generated  in  the  LSTF  had  smaller  heights  and  shorter  periods  as  compared  to  the 

288  SUPERTANK  waves.  Two  cases  with  distinctively  different  breaker  types,  one  spilling 

289  and  one  plunging,  were  examined. 
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290  The  spilling  wave  case  was  initiated  with  the  Dean  equilibrium  beach  profile  (Eq. 

291  8).  Because  of  the  smaller  wave  heights,  beach-profile  change  occurred  at  a  slower  rate. 

292  Similar  to  the  first  SUPERTANK  wave  run,  ST_10A,  a  subtle  bar  formed  over  the  initial 

293  monotonic  beach  profile  (Fig.  6  upper).  For  the  spilling  wave  case,  the  upper  limit  of 

294  beach  change  was  0.23  m  above  MWL,  as  identified  from  the  smaller  scale  plot  (Fig.  6 

295  lower).  The  profile  converges  on  the  seaward  slope  of  the  offshore  bar.  For  the  LSTF 

296  plunging  wave  case,  shoreline  advance  occurred  with  each  wave  run  along  with  sustained 

297  onshore  migration  of  the  bar  (Fig.  7).  The  accumulation  at  the  shoreline  was  subtle,  but 

298  can  be  identified  if  viewed  locally  (Fig.  7  lower).  The  upper  limit  of  beach-profile 

299  change  was  located  at  0.26  m  above  MWL,  with  the  lower  limit  identified  at  the  profile 

300  convergence  point  midway  on  the  seaward  slope  of  the  bar.  Overall,  the  trends  observed 

301  in  the  three-dimensional  LSTF  experiment  are  comparable  to  those  in  the  two- 

302  dimensional  SUPERTANK  experiment. 

303  Table  2  summarizes  the  upper  and  lower  limits  of  change  during  each  wave  run, 

304  including  the  breaking  wave  height.  In  summary,  for  the  30  SUPERTANK  wave  runs 

305  and  two  LSTF  wave  cases,  the  incident  breaking  wave  height  ranged  from  0.26  to  1.18  m 

306  (Table  2).  The  measured  upper  limit  of  profile  change,  including  the  scarped  dune  cases, 

307  ranged  from  0.23  to  0.70  m.  The  lower  limit  of  beach  change  ranged  from  0.50  to  1.61  m 

308  below  MWL.  Relationships  between  the  profile  change  and  wave  conditions  are 

309  discussed  in  the  following  sections. 


310 
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3 1 1  Cross-shore  Distribution  of  Wave  Height 

3 12  Wave-height  decay  is  representative  of  the  energy  dissipation  as  a  wave 

313  propagates  onshore.  Wave  decay  patterns  were  measured  by  the  closely  spaced  gauges 

3 14  for  both  the  SUPERTANK  and  the  LSTF  experiments.  Figure  8  shows  time-series  wave 

315  decay  patterns  measured  at  the  first  SUPERTANK  wave  run,  ST_10A.  As  discussed 

316  above,  considerable  beach  profile  change,  for  example  the  fonnation  of  an  offshore  bar, 

317  was  produced  during  this  wave  run  (Fig.  2).  The  substantial  morphology  change  also 

318  influenced  the  pattern  of  wave  decay.  The  point  of  steep  wave  decay  migrated  slightly  as 

319  the  beach  morphology  changed  from  the  initial  monotonic  profile  to  a  barred-beach 

320  profile.  This  point  was  defined  as  the  location  and  height  at  which  the  wave  breaks 

32 1  (Wang  et  al.,  2002).  For  ST_10A,  the  significant  breaking  wave  height  was  0.68  m. 

322  The  rate  of  wave-height  decay  tended  to  be  smaller  in  the  mid-surf  zone  (10  to  20  m)  as 

323  compared  to  the  breaker  zone  (20  to  25  m)  and  the  inner  surf  zone  (landward  of  10  m). 

324  The  offshore  wave  height  remained  largely  constant  until  reaching  the  breaker  line. 

325  The  wave  decay  pattern  for  the  longer  period  accretionary  wave  run,  ST_30A  (Fig. 

326  9),  was  considerably  different  than  the  steep  erosive  waves.  The  significant  breaking 

327  wave  height  was  0.36  m.  The  time-series  wave  pattern  remained  constant  for  each  wave 

328  run,  apparently  not  influenced  by  the  subtle  morphology  change  (Fig.  3).  The  bar  was 

329  formed  during  the  previous  wave  runs  with  greater  wave  heights.  Therefore,  instead  of 

330  breaking  over  the  bar,  shoaling  or  increase  in  height  of  the  long  period  wave  was 

331  measured  (at  around  30  m).  The  main  breaker  line  was  identified  at  around  15  m,  where 

332  a  sharp  drop  in  wave  height  was  observed. 
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333  For  the  dune  erosion  run,  ST  60A,  the  wave  height  remained  largely  constant 

334  offshore  (Fig.  10).  Significant  wave-height  decay  was  measured  over  the  offshore  bar  at 

335  approximately  30  m,  with  a  breaker  height  of  0.61  m.  A  noticeable  increase  in  wave 

336  height  was  measured  at  around  15  m  offshore,  likely  a  result  of  reflected  waves  off  the 

337  scarp,  followed  by  a  sharp  decrease  in  height  in  the  inner  surf  zone. 

338  The  cross-shore  distribution  of  wave  height  for  the  monochromatic  wave  run 

339  ST  IO  was  erratic  with  both  temporal  and  spatial  irregularity  (Fig.  1 1).  The  erratic  wave 

340  height  distribution  corresponds  to  the  irregular  beach-profile  change  observed  during  this 

341  wave  run  (Fig.  5).  The  breaking  wave  height  varied  considerably,  from  0.72  to  0.81  m, 

342  likely  caused  by  reflection  of  the  monochromatic  waves  from  the  beach  face.  The  wave- 

343  height  variation  in  the  offshore  region,  seaward  of  the  breaker  line  around  30  m,  was 

344  likely  related  to  oscillations  in  the  wave  tank. 

345  The  LSTF  experiments  were  designed  to  examine  the  effects  of  different  breaker 

346  types  on  sediment  transport  and  morphology  change.  Offshore  wave  heights  of  0.27  m 

347  were  generated  for  both  cases  (Fig.  12),  which  had  different  wave  periods.  However,  the 

348  cross-shore  distribution  of  wave  heights  was  considerably  different.  The  wave-height 

349  decay  at  the  breaker  line  was  much  greater  for  the  plunging  case  than  for  the  spilling  case, 

350  as  expected.  The  breaking  wave  height  was  0.26  m  and  0.27  m  for  the  spilling  and 

35 1  plunging  wave  runs,  respectively. 

352 

353  Wave  Runup 

354  The  extent  and  elevation  of  wave  runup  for  the  SUPERTANK  experiments  were 

355  measured  directly  by  the  closely  spaced  swash  gauges  (Kraus  and  Smith,  1994).  Figure 
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356  13  shows  the  cross-shore  distribution  of  time-averaged  water  level  and  wave  runup  for 

357  the  erosive  wave  run,  ST_10A.  The  swash  zone  water  level  was  measured  by  the  discrete 

358  swash  gauges,  as  discussed  above,  and  does  not  represent  time-averaged  water  level. 

359  Elevated  water  levels  were  measured  in  the  surf  zone.  As  expected,  the  mean  water  level 

360  in  the  offshore  area  remained  around  zero.  It  is  necessary  to  separate  the  elevation 

361  caused  by  wave  setup  and  swash  runup.  An  inflection  point  (labeled  with  an  arrow  and 

362  “ip”)  can  be  identified  from  the  cross-shore  distribution  curve  of  the  mean  water  level 

363  (Fig.  13).  The  inflection  point  also  tends  to  occur  around  the  still-water  shoreline  and  is 

364  regarded  here  as  the  distinction  between  wave  setup  and  swash  runup.  For  this  run,  the 

365  setup  measured  at  the  still-water  shoreline  was  0. 1  m,  which  is  about  17  percent  of  the 

366  total  wave  runup  of  0.6  m. 

367  For  the  accretionary  wave  run,  ST  30A,  the  inflection  point  in  the  mean-water 

368  level  also  occurs  around  the  still-water  shoreline  (Fig.  14).  The  average  setup  at  the 

369  shoreline  was  approximately  0.07  m,  also  about  17  percent  of  the  total  wave  runup  of  0.4 

370  m.  Total  wave  runup  was  significantly  limited  by  the  vertical  scarp  as  shown  in  the  dune 

371  erosion  run  of  ST  60A  (Fig.  15).  A  broad  setdown  was  measured  just  seaward  of  the 

372  main  breaker  line.  The  inflection  point  of  the  cross-shore  distribution  of  the  mean  water 

373  level  occurred  between  10  and  1 1  m,  before  reaching  the  still- water  shoreline  at  8  m.  The 

374  setup  measured  at  the  inflection  point  at  1 1  m  was  approximately  0.03  m.  The  wave 

375  setup  contributed  18  percent  of  the  total  wave  runup  of  0.17  m,  similar  to  the  above  two 

376  cases. 

377  The  cross-shore  distribution  of  time-averaged  mean  water  level  and  wave  runup 

378  for  ST_I0  (monochromatic  waves)  was  erratic  (Fig.  16).  As  opposed  to  the  irregular 
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379  wave  cases,  a  zero  mean-water  level  was  not  measured  at  a  considerable  number  of 

380  offshore  wave  gauges.  In  addition,  significant  variances  among  different  wave  runs  were 

381  also  measured.  The  total  wave  runup  varied  from  0. 16  to  0.35  m,  with  an  average  of  0.26 

382  m.  The  inflection  point  in  the  mean-water  level  distribution  occurs  around  the  still-water 

383  shoreline  at  5  m.  The  maximum  setup  at  the  still-water  shoreline  was  0.08  m,  which  is  3 1 

384  percent  of  the  total  wave  runup.  The  smaller  contribution  of  the  swash  runup  to  the  total 

385  wave  runup  can  be  attributed  to  the  lack  of  low-frequency  motion  in  the  monochromatic 

386  waves. 

387 

388  DISCUSSION 

389  Relationship  between  Wave  Runup,  Incident  Wave  Conditions,  and  Limit  of  Beach 

390  Profile  Change 

391 

392  The  measured  breaking  wave  height,  upper  limit  of  beach-profile  change,  and 

393  total  wave  runup  from  the  SUPERTANK  experiments  are  compared  in  Figure  17.  The 

394  thirty  runs  examined  are  divided  into  three  categories  describing  non-scarped  random 

395  wave  runs,  scarped  random  wave  runs,  and  monochromatic  wave  runs.  For  the  16  non- 

396  scarped  random  wave  runs,  except  the  three  runs  (10B20ER,  10E  270ER  and 

397  30D  40AR),  the  elevations  of  wave  runup  and  upper  limit  of  beach  change  roughly  equal 

398  the  significant  breaking  wave  height.  All  three  outliers  had  relatively  lower  measured 

399  swash  runup.  The  discrepancy  may  be  caused  by  the  performance  of  the  capacitance 

400  gauge.  The  partially  buried  capacitance  gauges  in  the  swash  zone  required  the  sand  to  be 

401  fully  saturated  (Kraus  and  Smith,  1994).  Both  10B  20ER  and  30D  40AR  are  initial 
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402  short-duration  runs,  during  which  the  sand  may  not  have  been  fully  saturated.  However 

403  this  does  not  explain  the  discrepancy  for  10E  270ER,  with  cause  unknown. 

404  For  the  scarped  random  wave  runs,  the  breaking  wave  height  was  much  greater 

405  than  the  elevation  of  wave  runup,  which  was  limited  by  the  vertical  scarp.  Because  the 

406  upper  limit  of  beach  change  was  identified  at  the  toe  of  the  scarp,  a  relationship  among 

407  the  breaker  height,  wave  runup,  and  beach-profde  change  is  not  expected  for  the  scarped 

408  random  wave  runs.  The  much  lower  wave  runup  by  monochromatic  waves  as  compared 

409  to  the  breaker  height  was  likely  caused  by  the  lack  of  low-frequency  motion.  No 

410  relationship  could  be  found  among  the  three  parameters  for  monochromatic  waves. 

411  Based  on  the  above  observations  from  the  SUPERTANK  data  with  breaking  wave 

412  heights  ranging  from  0.4  to  1.2  m  (Fig.  17),  a  simple  relationship  between  the  measured 

413  wave  runup  height  on  a  non-scarped  beach  and  the  breaker  height  is  found: 

414  *tw=l-0tfbi  (10) 

415 

416  The  average  ratio  of  i?tw  over  Hhs  for  the  16  non-scarped  wave  runs  was  0.93,  with  a 

417  standard  error  on  the  mean  of  0.05.  Excluding  the  three  questionable  measurements, 

418  10B_20ER,  10E_270ER  and  30D_40AR,  the  average  Rlw/f/\^  was  1.01,  with  a  standard 

419  error  of  0.02.  To  be  conservative  because  of  limited  data  coverage,  a  value  of  1.0  was 

420  assigned  in  Eq.  (10).  Caution  should  be  exercised  in  applying  Eq.  (10)  to  higher  waves 

42 1  than  the  range  examined  here. 

422  Comparisons  of  the  measured  wave  runup  with  the  various  existing  empirical 

423  formulas  (Eqs.  4,  6,  and  7)  and  Eq.  (10)  are  summarized  in  Figure  18  and  Table  3.  It  is 

424  recognized  that  Eq.  (4)  predicts  significant  runup  height,  whereas  Eqs.  (6)  and  (7)  predict 

425  2%  exceedence  of  runup.  The  measured  runup  Rtw  from  SUPERTANK  represents  a 
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426  maximum  value  of  total  wave  runup  (Eq.  10).  As  shown  in  Fig.  18,  previous  formulas 

427  did  not  reproduce  the  measured  wave  runup  at  SUPERTANK.  For  the  16  non-scarped 

428  wave  runs,  81%  of  the  predictions  from  Eq.  (10)  fall  within  15%  of  the  measured  wave 

429  runup.  In  contrast,  for  Eqs.  (4),  (6)  and  (7),  only  25%,  6%  and  13%  of  the  predictions, 

430  respectively,  fall  within  15%  of  the  measured  values.  Eqs.  (6)  and  (7)  under-predicted 

43 1  the  measured  wave  runup  significantly  for  the  erosional  cases,  but  over-predicted  runup 

432  for  the  accretionary  wave  cases.  The  discrepancy  is  caused  by  the  substantially  greater 

433  value  of  t,  for  the  gentle  long-period  accretionary  waves  than  for  the  steep  short-period 

434  erosional  waves  (Table  1).  Agreement  between  measured  and  predicted  values  was 

435  reduced  by  including  the  surf  similarity  parameter,  The  simpler  Eq.  (4)  developed  by 

436  Guza  and  Thornton  (1982)  based  only  on  the  offshore  wave  height,  more  accurately 

437  reproduced  the  measured  values  of  wave  runup  than  Eqs.  (6)  and  (7). 

438  Equation  (10)  was  applied  to  the  three-dimensional  LSTF  experiments  with  lower 

439  wave  heights  than  in  SUPERTANK.  Although  wave  runup  was  not  directly  measured  in 

440  the  LSTF  experiments,  it  is  assumed  here  that  the  total  runup  is  equal  to  the  upper  limit  of 

441  beach-profile  change,  a  reasonable  assumption  as  verified  by  the  SUPERTANK  data. 

442  For  the  spilling  wave  case,  taking  the  upper  limit  of  beach  change  at  0.23  m  as  the  value 

443  for  total  wave  runup,  the  breaking  wave  height  of  0.26  m  resulted  in  an  over-prediction  of 

444  13%.  For  the  plunging  wave  case,  the  upper  limit  of  beach  change  was  0.26  m,  which  is 

445  almost  equal  to  the  0.27  m  breaking  wave  height.  Therefore,  the  LSTF  data,  with  a  finer 

446  grain  size  (0.15  mm)  than  SUPERTANK  (0.22  mm),  support  the  new  predictive  equation 

447  (Eq.  10). 
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448  The  dependence  of  wave  runup  on  beach  slope  has  been  questioned  by  various 

449  studies.  Douglass  (1992)  re-analyzed  the  Holman  (1986)  dataset  underlying 

450  development  of  Eq.  (6)  and  stated  that  runup  and  beach-face  slope  are  not  well  correlated. 

45 1  Douglass  argued  that  beach  slope  is  a  dependent  variable  that  is  free  to  respond  to  the 

452  incident  waves  and  should  not  be  included  in  runup  prediction.  Sunamura  (1984)  and 

453  Kriebel,  Kraus,  and  Larson  (1991)  found  dependencies  of  beach  slope  on  wave  height 

454  and  period,  the  latter  reference  giving  a  predictive  formula  expressed  in  terms  of  the 

455  Dean  number  (Eq.  9).  Nielsen  and  Hanslow  (1991)  found  a  relationship  between  the 

456  surf  similarity  parameter  and  runup  on  steep  beaches.  However,  for  gentle  beaches  with 

457  slopes  less  than  0.1,  they  suggested  that  the  surf  similarity  parameter  was  not  related  to 

458  runup.  A  subsequent  study  by  Hanslow  and  Nielson  (1993)  conducted  on  dissipative 

459  beaches  of  Australia  found  that  maximum  setup  did  not  depend  on  beach  slope. 

460  In  practice,  beach  face  slope  is  a  difficult  parameter  to  define  and  determine. 

461  Except  for  Stockdon  et  al.  (2006),  a  clear  definition  of  beach  slope  is  not  given  in  most 

462  studies.  Stockdon  et  al.  defined  the  foreshore  beach  slope  as  the  average  slope  over  a 

463  region  of  two  times  the  standard  deviation  of  a  continuous  water-level  record.  In 

464  predictive  modeling  of  morphology  change,  relations  between  runup  and  foreshore  slope 

465  would  be  interdependent.  In  the  present  study,  the  slope  was  defined  over  the  portion  of 

466  the  beach  extending  roughly  1  m  landward  and  seaward  from  the  shoreline.  Substantially 

467  different  beach  slopes  can  be  obtained  by  imposing  different  definitions.  Inclusion  of  the 

468  beach  slope  in  predictive  relations  for  wave  runup  thus  adds  ambiguity  in  applying  such 

469  formulations. 
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470  Determining  offshore  wave  height  may  also  introduce  uncertainty.  In  most  field 

47 1  studies,  the  offshore  wave  height  is  taken  to  be  the  measurement  at  a  wave  gauge  in  the 

472  study  area.  Similarly,  in  this  study  it  is  taken  as  the  wave  height  measured  at  the  farthest 

473  offshore  gauge.  The  definition  of  an  offshore  wave  height  varies  between  studies,  in 

474  which  it  is  often  taken  at  whatever  depth  the  instrument  is  deployed  (Guza  and 

475  Thornton,  1981;  Guza  and  Thornton,  1982;  Holman,  1986).  In  addition,  under 

476  storm  conditions,  estimation  of  the  offshore  wave  height  may  not  be  straightforward 

477  (WANG  et  al.  2006). 

478 

479  A  Conceptual  Derivation  of  the  Proposed  Wave  Runup  Model 

480 

48 1  Swash  uprush  on  a  sloping  beach  is  often  approximated  using  a  ballistics 

482  approach  of  bore  propagation  (Baldock  and  Holmes,  1999;  Coco,  O’Hare,  and 

483  Huntley,  1999;  Larson,  Kubota,  and  Erikson,  2004;  Mase,  1988;  Suhayda,  1974). 

484  Most  derivations  are  based  on  the  bore  runup  model  of  Shen  and  Meyer  (1963)  and  the 

485  radiation  stress  formulation  of  Longuet-Higgens  and  Stewart  (1962).  In  the 

486  following  derivation,  a  similar  approach  is  adopted  to  examine  the  physics  foundation  of 

487  Eq.  (10). 

488  Assuming  a  normally  incident  wave  and  neglecting  longshore  variations,  the 

489  forces  acting  on  a  water  element  in  the  swash  zone  in  the  cross-shore  direction,  x,  (Fig. 

490  19)  can  be  balanced  as: 

r  hTt 

491  -/jgAxAyAz  sin  /? - /tAxAvL2  =  /tAxAvAz — -  (11) 

8  c )t 
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where,  p  =  density  of  water;  g  =  acceleration  due  to  gravity;  Ax,  Ay,  and  A z,  =  length, 
width,  and  height,  of  the  water  element,  respectively;  sin  (t  =  beach  slope; /=  friction 
coefficient;  and  Vx  =  velocity.  Eq.  (11)  can  be  reduced  to 

^  =  -gsin/i--^Vf  (12) 

dt  8Az 

Assuming  the  friction  force  is  negligible,  an  assumption  supported  by  experiments 
discussed  in  Komar  (1998),  Eq.  (12)  is  further  reduced  to 

dV 

— ^  =  -gsin/?  (13) 

at 

Influences  of  friction  and  infiltration  on  swash  motion  are  discussed  in  Puleo  and 
Holland  (2001).  Integrating  Eq.  (13)  with  respect  to  time,  yields 

K  =  K,  ~  &  s'n  P  (14) 

where,  V0  =  initial  velocity.  Integrating  Eq.  (14)  again  with  respect  to  time  gives  the 
swash  excursion,  x,  as  a  function  of  time,  t : 

x(t)  =  V0t-^-sm/3  (15) 

From  Eq.  (14),  the  maximum  uprush  occurs  at  a  time,  tmax,  when  the  velocity  becomes 
zero: 


max  •  n 

S’ sin  p 


with  a  corresponding  value  of  maximum  swash  excursion  of 

V2 


(16) 


*(hnax): 


2  g  sin  p 


(17) 


Assuming  a  small  and  planar  foreshore  slope,  then  tan  ft  ~  sin  /;,  and  the  elevation  of  the 
maximum  swash  uprush  Rsr  ma_x, becomes: 
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R 


:  x(tmJ  tan /3  - 


V2  V2 

- 2 - tan  B  ~  — 2- 

2gsin/?  2g 


(18) 


Eq.  (18)  suggests  that  the  maximum  elevation  of  swash  runup  is  not  a  function  of  beach 
slope  if  bottom  friction  forcing  is  neglected. 

The  initial  velocity  V0  can  be  approximated  by  the  velocity  of  the  wave,  C.  In 
shallow  water,  the  wave  velocity  is  limited  by  the  local  water  depth,  hi. 

V.-C-l*  (19) 

Assuming  a  linear  relationship  between  local  breaking  or  brreaking  wave  height,  Hm,  and 
the  water  depth,  hi 

Hhl  =  yh,  (20) 

where  y  =  the  breaker  index.  Eq.  (19)  then  becomes 

HUI 


K  ~cz  =  g- 


y 


(21) 


Substituting  Eq.  (21)  into  Eq.  (18) 


R  _  =  -2- 


K2  _  gHb,  _  Hh 


2  g  2gy  2y 


(22) 


Because  of  wave  height  to  depth  scaling  in  the  surf  zone,  it  is  reasonable  to  assume  that 
the  initial  velocity  V0  can  be  taken  at  the  main  breaker  line.  With  significant  breaker 
height  Hbs  Eq.  (22)  then  becomes 


4ma x=^  =  aHbs 


2  7 


(23) 


where  a  =  l/2y.  Eq.  (23)  indicates  a  linear  relationship  between  breaking  wave  height 
and  the  maximum  swash  runup,  supporting  the  findings  deduced  from  the  SUPERTANK 
experiment. 
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532  Kaminsky  and  Kraus  (1994)  examined  a  large  dataset  on  breaking  wave  criteria 

533  that  included  both  laboratory  and  field  measurements,  They  found  the  majority  of  y 

534  values  range  from  0.6  to  0.8,  which  yields  a  values  from  0.63  to  0.83.  Based  on  previous 

535  discussion  (Figs.  13  through  16),  swash  runup  constitutes  approximately  83%  of  the  total 

536  wave  runup.  Adding  the  17%  contribution  from  the  wave  setup,  the  total  wave  runup  Rtw 

537  is  roughly  equal  to  the  breaking  wave  height,  further  supporting  the  predictive  equation 

538  developed  from  the  SUPERTANK  dataset.  Thus,  the  empirical  model  of  total  wave 

539  runup  developed  based  on  the  SUPERTANK  data  is  supported  by  an  accepted  physical 

540  picture.  In  addition,  little  ambiguity  exists  in  the  straight-forward  parameterization  as 

541  given  in  Eq.  10. 

542 

543  CONCLUSIONS 

544  The  SUPERTANK  data  set  indicates  that  the  vertical  extent  of  wave  runup  above 

545  mean  water  level  on  a  non-scarped  beach  is  approximately  equal  to  the  significant 

546  breaking  wave  height.  A  simple  formula  for  predicting  the  total  wave  runup  Rt w  =1.0  Hbs 

547  was  developed  by  comparison  to  measurements  and  justified  by  a  derivation  based  on 

548  ballistic  theory  of  swash  motion.  This  formula  does  not  include  beach  slope,  which  is 

549  difficult  to  measure  and  is  itself  dependent  on  wave  properties.  The  new  model  was 

550  applied  to  the  three-dimensional  LSTF  experiments  and  accurately  reproduced  the 

55 1  measured  wave  runup.  Inclusion  of  the  slope-dependent  surf  similarity  parameter 

552  decreased  the  accuracy  of  the  calculated  wave  runup  as  compared  to  the  measured  values. 

553  An  exception  to  the  direct  relationship  between  breaking  wave  height  and  runup 

554  concerns  the  presence  of  dune  or  beach  scarping.  The  steep  scarp  substantially  limits  the 
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555  uprush  of  swash  motion,  resulting  in  a  much  reduced  maximum  level,  as  compared  with 

556  the  non-scarping  cases.  For  monochromatic  waves,  the  measured  wave  runup  was  much 

557  smaller  than  the  breaking  wave  height.  The  lack  of  low-frequency  modulation  limits  the 

558  wave  runup  for  monochromatic  waves. 

559  Based  on  the  SUPERTANK  and  LSTF  experiments,  the  upper  limit  of  beach- 

560  profile  change  was  found  to  be  approximately  equal  to  the  total  vertical  excursion  of 

561  wave  runup.  Therefore,  the  breaking  wave  height  can  be  used  to  provide  a  reliable 

562  estimate  of  the  limit  of  wave  runup  which,  in  turn,  can  serve  as  an  approximation  of  the 

563  landward  limit  of  beach  change:  UL  =  Rtw  =  Hbs .  Physical  situations  that  are  exceptions 

564  to  this  direct  relationship  are  beaches  with  beach  or  dune  scarps.  For  the  scarped  cases, 

565  the  upper  limit  of  beach  change  was  much  higher  than  the  total  swash  runup  and  was 

566  controlled  by  the  elevation  of  the  benn  or  dune. 

567 
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575 

576  APPENDIX  I  -  NOTATION 

577 

578  The  following  symbols  are  used  in  this  paper: 

579  A  shape  parameter  relating  to  grain  size  and  fall  velocity 

580  C  wave  velocity 

581  /  friction  coefficient 
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582 

g 

gravitational  acceleration 

583 

h 

still-water  depth 

584 

H 

wave  height 

585 

Hh 

breaking  wave  height 

586 

Hbl 

local  breaking  wave  height 

587 

Hhs 

significant  breaking  wave  height 

588 

Hb  h 

high  frequency  component  of  wave  height  at  the  breaker  line 

589 

Hb  i 

low  frequency  component  of  wave  height  at  the  breaker  line 

590 

h, 

local  water  depth 

591 

Ho 

significant  deepwater  wave  height 

592 

Hs\  h 

high  frequency  component  of  wave  height  at  the  shoreline  line 

593 

Hn 

low  frequency  component  of  wave  height  at  the  shoreline  line 

594 

Lb 

lower  limit  of  beach  change 

595 

L0 

deepwater  wavelength 

596 

n 

spectral  peakedness  parameter 

597 

N 

Dean  number 

598 

Rs 

significant  wave  runup 

599 

Rsr  max 

elevation  of  maximum  swash  uprush 

600 

Rtw 

total  wave  runup 

601 

Ri 

2%  exceedence  of  runup 

602 

T 

wave  period 

603 

t  max 

time  of  maximum  swash  excursion 

604 

T 

peak  spectral  wave  period 

605 

uL 

upper  limit  of  beach  change 

606 

Vo 

initial  velocity 

607 

vx 

velocity  of  a  water  particle  in  the  across  shore 

608 

w 

sediment  fall  velocity 

609 

X 

cross-shore  coordinate;  horizontal  distance  from  the  shoreline 

610 

tan  jS 

beach  slope 

611 

tan  jSf 

foreshore  beach  slope 

612 

y 

breaker  index 

613 

Ax 

length  of  a  water  particle 

614 

Av 

width  of  a  water  particle 

615 

A z 

height  of  a  water  particle 

616 

n 

wave  setup 

617 

nM 

wave  setup  under  monochromatic  waves 

618 

% 

wave  setup  at  the  shoreline 

619 

k 

surf-similarity  parameter 

620 

621 

p 

water  density 

622 
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734  FIGURES 

735 

736  Figure  1.  The  SUPERTANK  experiment  (upper)  and  LSTF  (lower)  during  wave  runs. 

737 

738  Figure  2.  The  first  SUPERTANK  wave  run,  ST_10A  erosional  case.  Substantial 

739  shoreline  erosion  occurred  on  the  initial  mono  tonic  profile  with  the  development  of 

740  an  offshore  bar.  The  horizontal  axis  “distance”  refers  to  the  SUPERTANK 

741  coordinate  system  and  is  not  directly  related  to  morphological  features. 

742 

743  Figure  3.  The  SUPERTANK  ST_30A  accretionary  wave  run.  There  was  subtle  beach 

744  face  accretion,  with  an  onshore  migration  of  the  offshore  bar  (upper).  The  accretion 

745  near  the  shoreline  is  identified  if  viewed  at  local  scale  (lower). 

746 

747  Figure  4.  The  SUPERTANK  ST_60A  dune  erosion  wave  run.  A  nearly  vertical  scarp 

748  developed  after  40  min  of  wave  action,  with  the  upper  limit  of  beach  change 

749  identified  at  the  toe  of  the  dune  scarp. 

750 

751  Figure  5.  The  SUPERTANK  ST  IO  accretionary  monochromatic  wave  run.  The 

752  resulting  beach-profile  under  monochromatic  waves  is  erratic  and  undulating. 

753 

754  Figure  6.  The  LSTF  spilling  wave  case.  Erosion  occurred  in  the  foreshore  and  inner  surf 

755  zone.  The  eroded  sediment  was  deposited  on  an  offshore  bar. 


756 
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757  Figure  7.  The  LSTF  plunging  wave  case.  Slight  foreshore  accretion  and  landward 

758  migration  of  the  offshore  bar  occurred  during  the  wave  run. 

759 

760  Figure  8.  Cross-shore  wave-height  distribution  measured  during  SUPERTANK  ST_10A. 

761 

762  Figure  9.  Cross-shore  wave-height  distribution  measured  during  SUPERTANK  ST_30A. 

763 

764  Figure  10.  Cross-shore  wave-height  distribution  measured  during  SUPERTANK 

765  ST60A. 

766 

767  Figure  11.  Cross-shore  wave-height  distribution  measured  during  SUPERTANK  ST  IO. 

768 

769  Figure  12.  Cross-shore  wave-height  distribution  for  the  LSTF  spilling  and  plunging 

770  wave  cases. 

771 

772  Figure  13.  Wave  runup  measured  during  SUPERTANK  ST10A.  The  arrow  with 

773  notation  “ip”  refers  to  the  inflection  point  between  the  wave  setup  and  runup. 

774 

775  Figure  14.  Wave  runup  measured  during  SUPERTANK  ST_30A.  The  arrow  with 

776  notation  “ip”  refers  to  the  inflection  point  between  the  wave  setup  and  runup. 

777 

778  Figure  15.  Wave  runup  measured  during  SUPERTANK  ST  60A.  The  arrow  with 

779  notation  “ip”  refers  to  the  inflection  point  between  the  wave  setup  and  runup. 
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780 

781  Figure  16.  Wave  runup  measured  during  SUPERTANK  ST  IO.  The  arrow  with  notation 

782  “ip”  refers  to  the  inflection  point  between  the  wave  setup  and  runup. 

783 

784  Figure  17.  Relationship  between  breaking  wave  height,  upper  limit  of  beach  profde 

785  change  and  wave  runup  for  the  thirty  SUPERTANK  cases  examined. 

786 

787  Figure  18.  Comparison  of  measured  and  predicted  wave  runup. 

788 

789  Figure  19.  Forces  acting  on  a  water  element  in  the  swash  zone. 
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Table  1.  Summary  of  Selected  Wave  Runs  and  Input  Wave  and  Beach  Conditions 
(Notation  is  explained  at  the  bottom  of  the  table). 


Wave  Run 

ID 

Ho 

m 

Tp 

s 

L0 

m 

n 

N 

Hbs 

m 

tan/3 

Hb_h 

m 

Hb_, 

m 

Hsi_h 

m 

Hs 

m 

SUPERTAI 

NK 

10A  60ER 

0.78 

3.0 

14.0 

20 

6.4 

0.68 

0.10 

0.42 

0.66 

0.15 

0.13 

0.24 

10A  130ER 

0.78 

3.0 

14.0 

20 

6.8 

0.68 

0.09 

0.38 

0.67 

0.15 

0.10 

0.23 

10A  270ER 

0.78 

3.0 

14.0 

20 

6.9 

0.68 

0.10 

0.42 

0.65 

0.16 

0.10 

0.24 

10B  20ER 

0.71 

3.0 

14.0 

3.3 

6.6 

0.65 

0.14 

0.58 

0.63 

0.17 

0.10 

0.23 

10B  60ER 

0.73 

3.0 

14.0 

3.3 

6.8 

0.67 

0.11 

0.44 

0.65 

0.17 

0.11 

0.24 

10B  130ER 

0.72 

3.0 

14.0 

3.3 

7.0 

0.69 

0.09 

0.36 

0.67 

0.18 

0.12 

0.25 

10E  130ER 

0.69 

4.5 

31.6 

20 

4.9 

0.72 

0.11 

0.69 

0.71 

0.15 

0.15 

0.16 

10E  200ER 

0.69 

4.5 

31.6 

20 

5.0 

0.74 

0.12 

0.77 

0.72 

0.15 

0.15 

0.18 

10E  270ER 

0.69 

4.5 

31.6 

20 

5.1 

0.76 

0.09 

0.58 

0.74 

0.15 

0.16 

0.20 

10F  110ER 

0.66 

4.5 

31.6 

3.3 

5.1 

0.75 

0.09 

0.58 

0.72 

0.18 

0.15 

0.26 

10F  130ER 

0.68 

4.5 

31.6 

3.3 

5.1 

0.76 

0.08 

0.48 

0.74 

0.18 

0.13 

0.21 

10F  170ER 

0.69 

4.5 

31.6 

3.3 

5.1 

0.76 

0.08 

0.50 

0.73 

0.20 

0.12 

0.24 

GO  60EM 

1.05 

3.0 

14.0 

M 

10.0 

1.18 

0.10 

0.43 

1.18 

0.01 

0.11 

0.03 

GO  140EM 

1.04 

3.0 

14.0 

M 

10.5 

1.04 

0.10 

0.41 

1.04 

0.04 

0.08 

0.10 

GO  210EM 

1.15 

3.0 

14.0 

M 

10.8 

1.07 

0.09 

0.39 

1.07 

0.04 

0.11 

0.02 

30 A  60AR 

0.34 

8.0 

99.9 

3.3 

1.6 

0.41 

0.14 

2.24 

0.40 

0.06 

0.24 

0.08 

30  A  130AR 

0.33 

8.0 

99.9 

3.3 

1.6 

0.39 

0.13 

2.09 

0.38 

0.06 

0.24 

0.09 

30 A  200AR 

0.34 

8.0 

99.9 

3.3 

1.6 

0.41 

0.13 

2.02 

0.40 

0.06 

0.25 

0.10 

30C  130AR 

0.31 

9.0 

126.4 

20 

1.4 

0.40 

0.13 

2.36 

0.40 

0.04 

0.18 

0.05 

30C  200AR 

0.31 

9.0 

126.4 

20 

1.4 

0.39 

0.15 

2.31 

0.38 

0.04 

0.19 

0.06 

30C  270AR 

0.31 

9.0 

126.4 

20 

1.4 

0.39 

0.15 

2.60 

0.38 

0.04 

0.20 

0.06 

30D  40AR 

0.37 

9.0 

126.4 

20 

1.4 

0.42 

0.13 

2.00 

0.42 

0.05 

0.17 

0.07 

10  80AM 

0.60 

8.0 

99.9 

M 

2.9 

0.76 

0.20 

2.78 

0.76 

0.01 

0.38 

0.03 

10  290AM 

0.63 

8.0 

99.9 

M 

3.1 

0.81 

0.17 

2.35 

0.81 

0.01 

0.34 

0.02 

10  590AM 

0.60 

8.0 

99.9 

M 

2.7 

0.72 

0.12 

1.64 

0.73 

0.01 

0.25 

0.03 

60 A  40DE 

0.69 

3.0 

14.0 

3.3 

6.2 

0.61 

0.12 

0.55 

0.58 

0.14 

0.16 

0.24 

60 A  60DE 

0.69 

3.0 

14.0 

3.3 

6.2 

0.61 

0.10 

0.46 

0.60 

0.14 

0.12 

0.24 

60B  20DE 

0.64 

4.5 

31.6 

3.3 

4.4 

0.66 

0.11 

0.74 

0.63 

0.15 

0.18 

0.24 

60B  40DE 

0.63 

4.5 

31.6 

3.3 

4.4 

0.66 

0.11 

0.76 

0.62 

0.16 

0.18 

0.25 

60B  60DE 

0.65 

4.5 

31.6 

3.3 

4.4 

0.66 

0.12 

0.79 

0.63 

0.17 

0.18 

0.30 

LSTF 

Spilling 

0.27 

1.5 

3.5 

3.3 

10.0 

0.26 

0.11 

0.41 

N/C 

N/C 

N/C 

N/C 

Plunging 

0.24 

3.0 

14.0 

3.3 

4.4 

0.27 

0.13 

0.96 

N/C 

N/C 

N/C 

N/C 

H0  =  offshore  wave  height;  Tp  =  peak  wave  period;  L0  =  offshore  wavelength;  n  =  spectral  peakedness;  N 
=  Dean  Number;  Hbs  -  significant  breaking  wave  height;  tan  /J  =  beach  slope  defined  as  the  slope  of  the 


section  approximately  1  m  landward  and  1  m  seaward  of  the  shoreline;  f  =  surf  similarity  parameter;  Hb_b 
=  incident  band  wave  height  at  the  breaker  line;  Hb_i  =  low  frequency  band  wave  height  at  the  breaker 
line;  Hsi_h  =  incident  band  wave  height  at  the  shoreline;  Hs i_i  =  low-frequency  band  wave  height  at  the 
shoreline;  M  =  monochromatic  wave;  N/C  =  Not  calculated. 


Table  2.  Summary  of  Beach  Change  and  Breaking  Wave  Height 


Wave  Run 

ID 

Hbs 

m 

uL 

m 

Ll 

m 

Scarp 

SUPERTANK 

10A  60ER 

0.68 

0.66 

1.29 

No 

10A  130ER 

0.68 

0.66 

1.29 

No 

10A  270ER 

0.68 

0.66 

1.29 

No 

10B  20ER 

0.65 

0.67 

1.35 

No 

10B  60ER 

0.67 

0.67 

1.35 

No 

10B  130ER 

0.69 

0.67 

1.35 

No 

10E  130ER 

0.72 

0.74 

1.52 

No 

10E  200ER 

0.74 

0.84 

1.52 

No 

10E  270ER 

0.76 

0.84 

1.52 

No 

10F  110ER 

0.75 

0.43 

1.52 

Yes 

10F  130ER 

0.76 

0.42 

1.52 

Yes 

10F  170ER 

0.76 

0.48 

1.52 

Yes 

GO  60EM 

1.18 

0.38 

1.61 

No 

GO  140EM 

1.04 

0.25 

1.61 

Yes 

GO  210EM 

1.07 

0.27 

1.61 

Yes 

30A  60AR 

0.41 

0.31 

1.36 

No 

30A  130AR 

0.39 

0.31 

1.36 

No 

30A  200AR 

0.41 

0.31 

1.36 

No 

30C  130AR 

0.40 

0.39 

1.01 

No 

30C  200AR 

0.39 

0.42 

1.01 

No 

30C  270AR 

0.39 

0.42 

1.01 

No 

30D  40AR 

0.42 

0.43 

0.65 

No 

10  80AM 

0.76 

0.46 

1.82 

No 

10  290AM 

0.81 

0.53 

1.82 

No 

10  590AM 

0.72 

0.53 

1.82 

Yes 

60A  40DE 

0.61 

0.28 

1.16 

Yes 

60A  60DE 

0.61 

0.28 

1.16 

Yes 

60B  20DE 

0.66 

0.38 

0.99 

Yes 

60B  40DE 

0.66 

0.38 

0.99 

Yes 

60B  60DE 

0.66 

0.38 

0.99 

Yes 

LSTF 

Spilling 

0.26 

0.23 

0.62 

No 

Plunging 

0.27 

0.26 

0.50 

No 

UL,  Ll  =  upper  and  lower  limit  of  beach  change,  respectively. 


Table  3.  Summary  of  Measured  and  Predicted  Wave  Runup. 


Wave  Run 

ID 

Hbs 

m 

Ktw 

m 

Eq  4 
m 

Eq  6 
m 

Eq  7 
m 

Eq  10 
m 

10A  60ER 

0.68 

0.60 

0.59 

0.43 

0.31 

0.68 

10A  130ER 

0.68 

0.60 

0.59 

0.40 

0.29 

0.68 

10A  270ER 

0.68 

0.70 

0.59 

0.43 

0.31 

0.68 

10B  20ER 

0.65 

0.33 

0.54 

0.49 

0.38 

0.65 

10B  60ER 

0.67 

0.70 

0.55 

0.42 

0.31 

0.67 

10B  130ER 

0.69 

0.64 

0.55 

0.36 

0.26 

0.69 

10E  130ER 

0.72 

0.77 

0.52 

0.53 

0.46 

0.72 

10E  200ER 

0.74 

0.78 

0.52 

0.58 

0.51 

0.74 

10E  270ER 

0.76 

0.45 

0.52 

0.47 

0.41 

0.76 

30A  60AR 

0.41 

0.41 

0.28 

0.70 

0.72 

0.41 

30A  130AR 

0.39 

0.43 

0.27 

0.64 

0.66 

0.39 

30A  200AR 

0.41 

0.42 

0.28 

0.64 

0.66 

0.41 

30C  130AR 

0.40 

0.40 

0.25 

0.67 

0.73 

0.40 

30C  200AR 

0.39 

0.42 

0.25 

0.73 

0.79 

0.39 

30C  270AR 

0.39 

0.42 

0.25 

0.73 

0.79 

0.39 

30D  40AR 

0.42 

0.23 

0.30 

0.69 

0.76 

0.42 

Bold  font  indicates  predicted  values  that  fall  within  15%  of  the 
measured  runup.  Rt w  =  total  measured  wave  runup. 


